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Abstract 

We propose in this paper a generic model of a non-standard aggregation mechanism for self- 
assembly processes of a class of materials involving the mediation of intermediates consisting of a 
polydisperse population of nano-sized particles. The model accounts for a long induction period in 
the process. The proposed mechanism also gives insight on future experiments aiming at a more 
comprehensive picture of the role of self-organization in self-assembly processes. 
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I. INTRODUCTION 



Nanophase materials are of key importance in a variety of fields from molecular and 
cellular biology to technological innovation. They also pose major theoretical challenges in 
view of their multiple scale dynamics, whereby microscopic level processes determine macro- 
scopic properties through the growth of initial fluctuations favoring in some way the specific 
material that will eventually emerge from synthesis. This switches on, in turn, a variety of 
self-organization phenomena associated with the coexistence of competing pathways. 

The processes we investigate in this work involve two steps. First, an initial species, 
N, is able to agglomerate into an intermediate form, X. We will primarily envisage the 
increase/decrease in the amount of X as being due to the attachment/detachment of units 
of N although more complex dynamics might be relevant in some cases. We expect that N 
and X will reach an equilibrium in a relatively short time. The second step of the process 
is the self-assembly of the final product, S, from X. By "self-assembly" we imagine that the 
units of X must not only join together to form S, but must undergo some sort of internal 
transition (such as a restructuring) which is quite slow. One of the main points of this 
work is to show that simply assuming a long characteristic time (i.e. a small rate constant) 
for this internal transition is not sufficient to give a long induction time. However, by 
including cooperativity - whereby the nascent S material templates or catalyzes the internal 
transition - arbitrarily long induction times can be achieved. A mathematical model of the 
contribution of cooperativity will be developed based on general physical considerations of 
how the process must proceed. 

The need to synthesize intermediate substances prior to the material of interest is also a 
well known phenomenon in chemistry, where it is typically also reflected by the appearance of 
an induction period^"-. Hierarchical self-assembly is also ubiquitous in biology. In particular, 
in a series of papers on the polymerization and crystallization of sickle cell Hemoglobin (HbS) 
-an issue of crucial importance for the pathophysiology of sickle cell anemia- it has been 
shown that metastable dense liquid clusters serve as precursors to the ordered nuclei of 
the HbS polymerization^"-. Induction times causing time delay have been observed and 
analyzed during this non-standard nucleation process of the HbS polymers. More generally, 
it is now known that homogeneous nucleation -most notably in nanophase materials such as 
proteins but apparently even in simple fluids - can involve a two step process starting with 
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the formation of more structured liquid-like droplets and/or clusters^"— from the solution 
followed by the development of crystalline order— li^. Here, the material in solution can be 
viewed as the population N of our scheme, the dense-liquid phase as the intermediate, X, 
and the solid as S. Cooperativity is again responsible for the induction period and time 
delays during the subsequent ordering corresponding to the processes N ^ X and X ^ S. 

A class of nanophase materials in which self-assembly appears to be accompanied by 
a rich variety of unexpected behaviors are synthetic zeolitesii^i^"— which find nowadays 
applications in a broad range of areas such as heterogeneous catalysis, petroleum refining, 
and microelectronicsi^ii^. In the context of zeolite formation mechanism, Silicalite-1 type 
zeolite^S in particular has been the subject of numerous studies^"—. 

The starting point of Silicalite-1 synthesis is a suspension of 3 nm silicate nanoparti- 
cles^S. The first period of ageing, at room or elevated temperature, is characterized by a 
rapid evolution of these nanoparticles. During this period, the average nanoparticle diameter 
increases to a slightly larger value of 5-6 nnt^^'^'^'^. On a structural level, the nanopar- 
ticle silicate network condenses, i.e. additional siloxane (Si-O-Si) bonds are formed^i^i^. 
Subsequently, still larger particles start growing, ultimately resulting in Silicalite-1 zeolite 
crystals. Convincing evidence shows that Silicalite-1 crystals grow by aggregation of the 5-6 
nm nanoparticles^. It was proposed that a fraction of the condensed 5-6 nm nanoparticles 
develops a zeolite crystalline framework. These intermediate nanoparticles can be identified 
as "nuclei" for zeolite crystal growth by nanoparticle aggregation^. The picture is then one 
of nanoparticles N giving rise to intermediates X which, in turn, form the (zeolite-crystalline) 
solid S, as illustrated in Fig. [1] Recent experimental liquid state ^^Si NMR and pH results^ 
of Silicalite-1 zeolite crystallization are presented and described in Fig|2l 

Fig|2]-ab shows the experimental distribution of silicon among different species in suspen- 
sion, determined with liquid state ^^Si NMR. There is a long induction period during which 
the amount of nanoparticles varies slowly, followed by a rather abrupt transition towards an 
equilibrium state dominated by zeolite crystalline material (solid). Because nanoparticles N 
and intermediates X have similar size NMR cannot distinguish between these two species. 
The rapid condensation of nanoparticle framework in the initial stage, resulting in the for- 
mation of intermediates, is evident from the silicon connectivity distribution (Fig [2]-d) and 
the corresponding increase of the suspension pH (Fig[2]-c). Upon transition to crystalline 
material, the pH shows a second increase. 
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FIG. 1: A schematic demonstration of multistep versus standard synthetic pathways in a two-order 
parameter space accounting for the presence of partiahy structured intermediates. 

Notice that during the later stages of crystaUization, other mechanisms, such as Ost- 
wald ripening, acting on the aheady formed aggregates are hkely to take over as dominant 
crystaUization mechanism. 

In view of the foregoing we argue that a timescale separation and the appearance of 
long quasi- stationary plateaus is likely to be a generic property of large classes of synthetic 
processes. One of the main goals of this article is the theoretical explanation of the mecha- 
nisms underlying these phenomena. We establish a generic theoretical model which captures 
the main features of self-organization in the process of multi-step nucleation and growth. 
Specifically, we show how the combination of intermediate steps in self-assembly, together 
with first-order nucleation and cooperativity, in the form of growth controled by the avail- 
able surface area of the product, can explain the long induction times often observed in 
these processes. For the particular case of zeolites, our model is intended to sort out the 
basic microscopic mechanisms that control the nucleation and growth of the final crystalline 
product. This can have some practical repercussions in not only eliminating superfluous 
conjectures that could be entertained in the absence of theoretical input, but also in further 
designing targeted and better controlled experiments. 

The article is organized as follows: In Section 2, we outline the proposed mechanisms for 
self-organization incorporated in our mathematical models. In Section 3 we develop a simple 
phenomenological model based on the assumption of cooperative behavior which reproduces 
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FIG. 2: (Color on line) Liquid state ^^Si NMR and pH results from Ref. [36(]. (ab) Evolution of 
nanoparticles (N and similarly sized X particles), oligomers and NMR- undetected phase (denoted 
as "solid"). The undetected phase consists of large particles, i.e. crystalline solid (S*) and possibly 
a fraction of intermediates (X particles that are considerably larger than nanoparticles), that fall 
outside the NMR detection limit. The oligomers are an additional group of smaller silicate species, 
but their presence is assumed to be of only secondary importance for the crystal growth. Panel-a: 
direct NMR results obtained with a short delay time of 7 sec between radiofrequency pulses in 
the NMR experiment. The process shows two-step evolution. The first step (< Ihr) indicates the 
formation of intermediate species. Subsequently, there was a long induction period during which the 
amount of nanoparticles varied slowly, followed by a rather abrupt transition towards an equilibrium 
state dominated by zeolite crystalline material. It is noted that in the original publication, these 
NMR data were corrected to account for a longer NMR delay time of 92 s (Panel-b), based on 
measurement of the unheated sample^. The correction resulted in quantitatively more reliable 
populations but the initial step was suppressed. The two-step evolution was also evident from 
measurements of the suspension pH (panel-c) and the connectivity of the nanoparticle framework 
(panel-d). The connectivity is expressed as fractions, where stands for a silicon atom 
connected to n = 1 — 4 neighboring silicons via siloxane bonds. The disappearance of the 
fraction might indicate a depletion of the more structured X-particles. 
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much of the quahtative behavior summarized above. A more microscopic and mechanistic 
model is presented in Section 4 and shown to describe the generic quahtative features of 
recent observations and lead to a prediction regarding the distribution of cluster sizes. The 
main conclusions are summarized in Section 5. 

II. THE MECHANISMS FOR SELF-ORGANIZATION AND THEIR MODELING 

In this section we outline the construction of mathematical models capable of quantifying 
the mechanisms summarized in the Introduction and, in particular, the presence of multiple 
time scales and the existence of a long induction period. The models to be developed are 
articulated around three basic steps which account for the self-organizing behavior. 

(a) Equilibration between nanoparticles and intermediates. This part of the dynamics 

will depend on the details of specific systems. Here, we will model it as a globally 
autocatalytic process whereby nanoparticles and intermediates interact to form new 
intermediates while depleting the nanoparticle population. The rationale behind this 
assumption is that in the particular case of zeolites, the intermediate is not a single 
species but, rather, it is believed to consist of a population of structures of various 
sizes2ii2£. We assume that the dominant mechanism for conversion of nanoparticles 
into intermediates is growth of smaller intermediates via aggregation of nanoparticles, 
i.e., via an autocatalytic mechanism. It would be possible to include nucleation of new 
intermediates, but it seems likely that growth will be the dominant process so that, 
for simplicity, nucleation can be neglected. It is assumed that at the beginning, a 
small population of intermediates is already present due, e.g., to thermal fluctuations. 
Obviously, for other systems, such as proteins, first-order reactions characteristic of 
nucleation would be more appropriate. 

(b) Initial formation of final product. At some point, a stable cluster of S must be formed 

which will then proceed to grow and we assume that the long plateau referred to earlier 
is a measure of the time required for this to occur. In a simple fluid this would be 
recognized as a nucleation phenomena^ which is typically modeled as a probabilistic 
process occurring with a probability per unit time proportional to the supersaturation. 
In the case of more complex materials, it may be that the intermediates must pass 
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through multiple structural phases or that some other type of transformation of the 
intermediates is necessary before stable clusters of final product can form. It might 
be more appropriate in this case to speak of self-assembly rather than nucleation^. 
Nevertheless, we will assume that this complex microscopic process can be crudely 
described in the same terms as nucleation in a simple fluid. 

(c) Growth and Cooperativity. Once the product begins to form, the speed of the process 
increases dramatically. This is a signal of cooperative behavior whereby the presence 
of some product accelerates the formation of new product. One possibility is that 
this is simply due to the growth of the clusters once they are formed. In simple fluids, 
secondary nucleation - the nucleation of new crystals on the surfaces of existing crystals 
- also occurs and can speed the formation of crystals^E. In more complex materials, 
if the intermediates must pass through different phases, then a growing cluster could 
also catalyze that process thereby accelerating cluster formation. We will not try to 
distinguish these different mechanisms but will simply account for an acceleration of 
the rate of cluster formation due to the presence of existing product. 

We will discuss the implementation of this phenomenological framework in two stages. 
First, we translate it into a crude phenomenological model which serves to illustrate the 
general ideas. Then, we describe a more microscopic, mechanistic approach that leads 
to similar results and that has the advantage of being less ad hoc and of making further 
predictions. 

III. PHENOMENOLOGICAL MODEL 

We consider a closed reactor initially containing nanoparticles, N, present in abundance, 
along with an intermediate (more ordered and larger) species, X, present in small amounts. 
Rather than forming the product, S, directly, the N particles first accumulate to form the in- 
termediate X. In fact, the single interemediate in the model represents a spectrum of species. 
These species undergo reactions which convert from one type of intermediate to another, 
but this will not be modeled here. Instead, we only track the net mass of intermediates 
which increases by consuming nanoparticles. In certain cases, such as zeolite crystallization, 
it is believed that the intermediates are the only species possessing the necessary structure 
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that allows them to bind together to form solid^. For this reason, we restrict our analysis 
to the sequence N ^ X ^ S and do not consider possibilities such as direct conversion of 
nanoparticules into solid which may be of relevance to other processes. Such alternatives 
could easily be included within the framework described below. The N to X transition is 
taken to be an autocatalytic process of order higher than one, with a rate ui depending on 
the abundance of both N and X. Furthermore, in agreement with available data it is stip- 
ulated that the X particles accumulate to give rise to S. The X to S transition occurs at a 
rate, 1^2 , that describes two processes. First, as discussed above, the initially slow formation 
of stable clusters is modeled as a probabilistic event occurring at a rate that depends on the 
supersaturation. Then, as the final product is formed, the rate increases due to cooperativ- 
ity. Thus, the total rate of formation of the product will depend on the concentration of 
intermediates and of the product itself. To avoid the presence of stoichiometric coefficients 
we choose to work with the mass fractions n, x and s of N, X and S particles respectively 
since by mass conservation the rate of loss of n in the N to X transition will then necessarily 
be equal to the rate of growth of x, etc. We obtain in this way the following three coupled 
differential equations, 

dn 

- = -u,{n^x) 

— = ui{n,x) - 1/2(3;, s) 

^ = z/2(a;,s), (3.1) 
the corresponding kinetic scheme being 

N ^ X ( rate z/i) 

X — > S ( rate z/2) (3.2) 

Notice the conservation condition 

n + X + s = constant (3.3) 

This relation does not involve the mass fraction of oligomers as they are believed to play a 
secondary role in the process and are therefore neglected^^. 

To illustrate this behavior, we choose to model the rate ui as a Verhulst type growth 
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process^i^. 



i'i{n,x) = kinx — k2x'^ (3.4) 



whereby nanoparticles and intermediates react to produce more intermediate species. Notice 
that balance is achieved when n = so that the ratio of rate constants can be fixed by the 
quasi-equilibrium level of the mass fractions which occur at short times. The value of, say ki 
is then determined by the time the system takes to reach this quasi-stationary state. Thus, 
these parameters are completely determined by the short-time behavior of the system. 

The rate 1^2 is broken into two parts: 2^2(2;, s) = h'2i{x) + 1^22(2;, s). The first part depends 
only on the amount of X present and represents the nucleation/self- assembly process. For 
this, we take the simplest reasonable form, 1^21(2;) = k^i^x — Xq), where Xq is the equilibrium 
concentration of X so that x — xo represents the supersaturation. The second part of the rate 
represents the effect of cooperativity. We have not found it possible to model the behavior 
observed in the zeolite data using higher order terms involving only the concentration of 
intermediates. Instead, we have found it necessary to model this as an enhancement of 
the nucleation rate that depends on the amount of product already formed, 2^22(2;, s) = 
k^as'^i^x — Xo). In the following we always assume that xq = which implies that once the 
product is formed from X, it will not dissolve again. 

Putting these pieces together, the phenomenological model is given by 

— = —kinx + k2X^ (3.5) 

— ^ = kinx — k2X^ — /^sxfl + as^) 
at 

f- = ksxil + as^") (3.6) 

Long induction times are expected to arise when the cooperative part of the process is slow, 
which is to say that k^ << ki,k2- In this case, the model exhibits dynamics on two different 
time scales: at short times, the cooperative part can be neglected with the result that 
c = n + X is conserved (since virtually no product forms) and the model can be integrated 
to give n{t) = c — x{t) and 

^' kic - {ki + k2)x{0) {e->'^^' - 1) ^ ' 

For kict >> 1, this leads to an equilibrium with x = -j^^^- Thus, assuming that c is on 
the order of 1, we find equilibration on a time scale determined by ki and with relative 
fractions of material in the form of N and X determined by k2. This short-time equilibration 
is expected to describe the beginning of the self-assembly process when materials are first 
mixed, heat is first applied, etc. 



At longer times, nucleation and cooperativity contribute causing the system to evolve 
away from the short-time steady state. There is then a new, long-time steady state in which 
X — > 0. At intermediate times, the model must be solved numerically. Figure [3] shows the 
result of such a solution starting with initial conditions x(0) = 10~^ and ?7,(0) = 1 — a;(0), 
corresponding to an initial state in which almost all the mass is in the form of nanoparticles. 
We take k2 = 9ki so that the short-time equilibration leads to n = 0.9 and x = 0.1. To 
illustrate the presence of a long induction time, the value = 0.002A;i was used so that 
if 1/ki is on the order of a minute, then significant self-assembly of the final product only 
occurs after a time on the order of l/k^ or about 8 hours. The remaining constants, a and 
fi, primarily control the speed with which the system makes the transition from its short- 
time steady state to the final long-time steady state. These values were chosen somewhat 
arbitrarily to be k^a = ki and /i = 1.75. The figure shows that after a rapid equilibration, 
the system experiences a long period of slow nucleation that is followed by a rapid transition 
to the final steady state. The final transition begins to occur around kit ~ 250 correspond 
to t ~ 4hr for ki = 1/minute. Qualitatively, this behavior is similar to that observed in the 
zeolite experiments of Aerts et al^^ (compare to Fig. [2]). 

Finally, we note some specific features of this model that are motivated by the experi- 
mental data shown in Fig. |2l First, there is no evidence of nucleation of intermediates from 
nanoparticles. This is in part based on the fact that the nanoparticle concentration reaches 
a finite plateau even though the intermediates are depleted. Nucleation may a priori occur 
on much longer timescales than that of the experiment in which case the observed plateau 
is only quasi-stationary. This would then suggest adding a term —kiU to the equation for 
dn/dt and a corresponding term in the equation for dx/dt. However, the rate constant 
would be very small thus justifying our neglect of this term on these timescales. Second, 
the conversion of intermediates to solid is irreversible. This is motivated by the fact that 
dissolution of zeolite crystals in alkaline medium does not lead to the extraction of silicate 
species of the size of nanoparticles or larger— (which also explains the depletion of X (see 
fig. 2d)). This leaves open the possibility of an equilibration between the solid and the 
oligomers, but, as can be seen from Fig. |2l such a reaction, if it exists, is of secondary 
importance and must take place on longer time scales than that of the formation of solid 
from nanoparticles. Besides, the fact that the crystals eventually grow to macroscopic size 
suggests that the formation of solid is, for all practical purposes, irreversible. 
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FIG. 3: (Color on line) Solution of the phenomenological model, Eq. (l3.5p . with the rate constants 
given in the text. The full line is the mass fraction of the nanoparticles, the dotted line is that of 
the intermediates, x{t), and the broken line is that of the final product, s{t) 

IV. A MORE MICROSCOPIC APPROACH 

Having shown that the combination of nucleation/self-assembly and cooperativity appears 
to give a reasonable description of the mechanisms of self-assembly, we now describe a more 
concrete realization of this picture. It is expected that assumptions made here are both 
more physically plausible than those made above and more intuitively appealing. They are: 

1. The N ^ X reaction is the same as assumed in the previous section. 

2. Formation of new stable clusters of final product occurs at a rate proportional to the 
supersaturation. 

3. The initial size of a stable cluster is always the same. 

4. Once formed, clusters grow at a rate proportional to their surface area and to the 
supersaturation. 

5. All clusters are geometrically identical. 

6. The formation rate is enhanced by a factor proportional to the surface area of the 
existing clusters and to the supersaturation. 

Assumptions 1 and 2 are the same as made previously. Assumption 3 is perhaps the most 
problematic: in a simple fluid it might be justified on the grounds that growing clusters are 
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presumed to start out as a stable cluster which has a fixed size (at least for a given super- 
saturation). This may be true here. If not, it can be imagined that there is a distribution 
of initial cluster sizes and that we are describing an average over an ensemble of such clus- 
ters. Assumption 4 explicitly introduces cluster growth which was neglected in the previous 
section (or, at least, it was crudely lumped together with cooperativity). Assumption 5 is 
merely a convenience so that we can be specific concerning geometric factors such as the 
ratio of volume to surface area of a cluster. Assumption 6 seems the simplest, physically 
plausible way to introduce the idea of cooperativity. 

These models involve the area. A, and mass or, equivalently, volume, V, of growing 
clusters. For spherical clusters, we of course expect that A R'^,V ^ where R is the 
radius. In fact, these scaling relations should hold for three dimensional growth, with R some 
measure of linear size, regardless of the precise shape of the clusters. However, other modes 
are possible and are of relevance for growth of some protein molecules. For example, suppose 
that the cluster has the shape of a cylinder with crossectional area vrr^ and length L. Quasi- 
one-dimensional growth would occur if the cylinder lengthens with constant cross section. 
This would involve attachment on the fiat faces so that the active growth area would be 
27rr^, a constant, while the volume would be 7rr^L(t). Quasi- two-dimensional growth would 
occur if attachment happened on the curved part of the cylinder so that the crossectional 
area increases with time while the length remains constant. Then, the active growth area 
would be 27rr (t) L and the volume vrr^ (t) L. All of these cases can be summarized by the 
general scaling A = GaR^~^,V = GyR^ where R is the appropriate linear length scale 
and Ga and Gy are constant geometric factors. We will use this as our general model with 
the additional provision that D is an integer. We therefore rule out fractal growth patterns 
which could a priori be treated within our general framework but which do not admit of the 
mathematical simplifications we utilize below. 

To translate this into a mathematical model for growth, some new quantities must be 
defined. Let N'{t) be the number of clusters growing at time t and let m{t,t') be the mass, 
at time t, of a cluster nucleated at time t'. Note that by assumption 3, m(t, t) = mo for some 
constant mg. Finally, the total mass of nanoparticles, intermediates, solid and oligomers is 
denoted M, and is a constant throughout the experiment. The total mass of solid at time t, 
Ms{t), is simply given by a sum over the masses of all clusters present at time t. A cluster 
at time t has mass m(t, t') if it was nucleated at time t'. The number nucleated at time t' is 
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simply dt' so 



sit) = f m{t,t')^^^dt' {U 



dt' 

where we assume that there is no final product present at the beginning of the experiment, 
s(0) = 0. Thus, based on assumption 4, the mass of a cluster grows as 

m(t,t)=mo (4.9) 

for some constant, kg. Finally, the rate of appearance of new clusters is the sum of a 
nucleation-like term proportional to the supersaturation (assumption 2) and a cooperative 
term proportional to the total area of the clusters (assumption 6), 

f^ = M*„. + *..£^(M')^A' (4.10) 

Note that we include a factor of the total mass in the nucleation term. This is because the 
number of clusters is an extensive quantity - as is the total area calculated in the second 
term. The number of nuclei generated per unit time must also be extensive, hence the factor 
of the mass. 

At this point, the model seems much more complicated than the phenomenological model 
given above. However, after some simple manipulations (see Appendix), it can be described 
by a set of three first order differential equations which uses the same number of parameters 
as in the phenomenological model: 

— = —kinx + k2X^ (4-11) 

dx . ,2 d'siu) 

— kiux — k2X — '-fx- 



dt du 
du 

-dt=^'' 

Here, uif) is an auxiliary function with the initial condition m(0) = 0. The mass fraction of 
S is given by s(t) = s (n(t)) where 

s (u) = f 1 + — ^ V a, exp iXiu) - D-u (4.12) 
\ a J ^-^ a 

The constants Aj are the roots of the £)-th order equation 

"^^^-^fj^rf^"-' (4-13) 
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while the coefficients Oj are 

^ "^n(A.-A,)'''-' 

The case of general values of D is discussed in the Appendix. All of this is written in terms 
of the two dimensionless constants 

a = — , P = — (4.15) 

and the rate constant 

where mg = pGyRo- Thus, despite the fact that the physical model involves 3 parameters, 
kg,kn and ka, as well as the initial mass of a stable cluster, mg, and the total mass M and 
the density p, the data can be fit with only three parameters, Q!,/3 and 7, which is no more 
than the phenomenological model. 

The present model is now written in terms of rate equations which are exactly equivalent 
to eqns (14.8114. lOl) . Still, it does not have the same form as the phenomenological model 
proposed in the previous Section. It is, however, possible to make a correspondence between 
them. Since both the mass fraction of S, s{t), and the auxiliary variable u{t) increase mono- 
tonically with time, and since they are related by a simple algebraic equation, eqn (14.121) . 
we can imagine inverting that equation to get u as a function of s. An analytic inversion 
is not possible but for small u and s, this can be done perturbatively giving m as a power 
series in s. (See the Appendix for details.) Substitution into eqn (14. lip and specializing to 
D = 3 gives 

— = —kinx + k2x'^ (4-17) 
at 

dx , ,9 ^ f ^ a + 3 a + 32 



ds ( a + 3 a + 32 



which is very similar to the phenomenological model except that the single term is replaced 
by an infinite series in s. Thus the adjustment of p is seen as a way of accounting for the 
fact that the infinite sum is being approximated by a single term. 
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FIG. 4: (Color on line) Solution of the microscopic model, Eg. (14. Ill) , (lines) with the rate constants 
given in the text. The full line is the mass fraction of the nanoparticles, the dotted line is that of 
the intermediates , x{t), and the dashed line is the mass fraction of the solid, s{t). 

Figure m shows the result of a numerical solution of the model, Eg. (14. lip . The constants 
ki and k2 have the same values as used previously. The other constants, a = 8750, /3 = 
1.67 X 10"'^ and 7 = 6.55 x 10~^/ci, were chosen to give approximately the same point of 
crossing of the curves and similar final states as seen in the phenomenological model. In 
this case, the plateau resulting from the long induction time is even more pronounced and 
the transition to the steady state is even sharper than in the phenomenological model. 

To test whether or not cooperativity is really required, we have attempted to model similar 
behavior with ka = 0. Typical results are shown in Fig. [5l In this case, there are only two 
parameters and we find that it is not possible to simultaneously capture the long plateau 
up to kit ~ 200 and to capture the crossing of the curves at about kit ~ 310. Furthermore, 
in contrast to both the data and the model with cooperativity, the mass fractions do not 
show plateaus at long times. As can be seen in the figure, a decrease in /3 by a factor of 
10~^, with 7 adjusted to fix the point at which the curves cross, has only a small effect on 
the quality of the fit to the data. Further decrease in /3 has little effect on the length of the 
plateau. We conclude that the cooperativity is required to explain the presence of a long 
induction time. 

One advantage of the microscopic approach is that, as discussed in the Appendix, the 
size distribution of the clusters is easily calculated. The number of clusters with radius R 
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FIG. 5: (Color on line) Comparison of the model without cooperativity and the data. The solid 
and dotted lines are the nanoparticle mass fractions and the remaining curves are the solid mass 
fraction. The solid and dashed lines were calculated using (3 = 5.0 and 7 = 0.0.0034A;i; the dotted 
and dash-dotted lines correspond to /5 = 5 x 10~^ and 7 = 2.66A;i. Decreasing /5 gives a slightly 
longer plateau but never as long and sharp as that possible when cooperativity is present. 

or smaller at time t, M {R,t), is given by 

^0, if R<Ro 

iiRo<R<Ro {u{t) + 1) 

^ATKt)), if i? > Kt) + 1) i?o 



U{R,t) 



The function J\f {u) is related to the total number of clusters by M {t) = ^J\f (u{t)) and is 
closely related to s'(m). Its explicit form is 



D 



u 



Qj exp (Xju) . 



(4.18) 



The factor i?o {u (t) + 1) occurring in the size distribution is the radius of the largest possible 
cluster at time t and the prefactor — is the maximum number of clusters that could be 

^ mo 

created from the given amount of material. 



V. CONCLUSION 



We have shown that a self-organization process involving self-assembly/nucleation, 
growth and cooperativity mediated by intermediates consisting of a polydisperse popula- 
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tion of nanosized particles with varying structure accounts for the kinetics of self-assembly 
observed in a class of nano-sized materials. In particular, the models reproduce a long quasi- 
steady state and the long waiting times observed as well as the final apparent equilibrium 
between nanoparticles and clusters of final product. Note that the reason for the termination 
of the self-assembly process in our model is due to the depletion of the particles, X. At some 
point, the N ^ X transitions can not compensate for the fast incorporation of X into the 
final phase. As new particles can only be formed when some X is already present, the whole 
process stops when X is fully consumed. The model that we have introduced may account 
for a broad class of aggregation processes. There is evidence that the complex mechanism of 
zeolite synthesis, at least some features of it, can be described by this model. The two-step 
behavior of Fig. [2]-a resembles the evolution of article populations as depicted in Figs. E] and 
m especially if we assume that some part of the (larger) intermediates X belong to the un- 
detected (sohd) phase. Moreover, the extinction of structured nanoparticles/intermediates 
which announces the end of the growth process, seems confirmed by experiment (compare 
in Fig. |2]-d with X{t) in Figs. [3] and HI). 

The basic physical mechanisms were modeled at two levels. First, a purely phenomeno- 
logical model was presented which captures the basic ideas in their simplest form. Then, a 
more detailed and mechanistic picture was given. This involved new concepts, such as the 
distribution of cluster sizes, but in the end was reduced to a system of rate equations similar 
to the first model. In fact, it was noted that the phenomeno logical model could be viewed 
as a systematic approximation to the more detailed model. 

It was shown that the concept of cooperativity plays a key role in explaining the obser- 
vations. When cooperativity is removed from the model, it is not possible to simultaneously 
reproduce the long waiting time at the beginning of the process and the rapid growth that 
occurs after long times. 

One point that deserves emphasis is the limited nature of the assumptions underlying 
our model. It is based on analogies to the familiar ideas of crystallization^^ but we do 
not claim that the processes are necessarily the same. The nucleation event that initializes 
crystallization in a simple fluid is probably a much simpler process than the self-assembly 
in more complex materials such as zeolites. However, it makes sense that in both cases the 
rate at which the process happens must surely, in a first approximation, be proportional to 
the supersaturation. Similarly, the mechanism of cooperativity in a simple fluid is simply 
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the accelerating growth of crystalhtes as they become larger. In some cases, secondary 
nucleation on growing crystals is also a factor. In zeolites, however, it may be a true catalysis 
of the transformation of intermediaries from one form to another. Whatever the nanoscale 
mechanism, it is plausible that in a first approximation the rate must be proportional to the 
surface area of existing clusters and to the supersaturation. Our results strongly suggest that 
some form of cooperativity, whose specifics will depend on the system at hand, is necessary 
when long induction times precede rapid transformation from one state to another. 

It would be interesting to compare the results of our model with other experimental datai^^ 
and theoretical models that have been proposed^ for zeolites. Moreover, it is hoped that 
this model will provide a framework for pursuing further experiments. In particular, the 
model makes a definite prediction concerning the distribution of cluster sizes and it would 
seem that this could be checked via experiment. 
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Appendix A: Elaboration of the model of Section IIVI 
1. Derivation of eqns (|4.im4.16]) 

The model for the formation of solid can be summarized as 




dm{t, t') 



kgxA{t, t') 



dt 




(A.l) 
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The mass and surface area of a growing cluster will be related in a system-specific way. For 
example, if the cluster is spherical, then the area will scale as the square of the radius and 
the mass as its cube. On the other hand, if the growth is primarily one-dimension, e.g. a 
cylinder which grows by becoming longer, then the mass and area will both increase in time 
as the first power of the length. We will assume in general that there is some length scale, 
t') characterizing the size of clusters at time t' which have nucleated at time t and that 
the area scales as A{t,t') = G AR^~^{t,t'), where Ga is a geometric factor and D is the 
dimensionality of the growth process so that the mass will scale as m(t,t') = pGvL^(t,t'). 
Then, the equations for the mass and number become 
dR (t, t') kgGA 



dt DpG, 



-X, R{t,t) = Ro (A.2) 



MknX + GAkaX R^{t, t')^^^^j^dt' 



Note that the initial radius is related to the initial mass of a cluster by where ttt-q = pGyR^ 
Next, introduce the variable u{t) = !MM^3^_ From 

i?(t,0 = i?o + ^ fx{t")dt" 

J-^P'^v Jt' 

we derive the following expression 

R{t, t') = R{t, 0) - R{t', 0) + Ro = {u{t) - u{t') + 1) Ro 



D 



From 



and 



du _ 1 dR{t, 0) _ kgGA 
dt Ro dt DpGy 



x(t) = jx(t) 



dx ds ds du 

we have already derived eqns |4.11] and |4.16] . However, we still need to prove eqns (14.121 - 
MB- 

Notice that the radius increases monotonically with time. We can therefore replace time 
by u{t) and write A/"(t) = ■^J\f{u(t)) for some function A/", (the prefactor being introduced 
to simplify later expressions). Using ^ = — the equation for the number of clusters 



can be rearranged to give 



^ - P + a {u-u' + 1) \ ' du' (A.3) 



du Jo du' 
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with 



DGykapRo DkaTriQ moDGvpRokn 

OL = ; = ; , p = - 



kn 



G.k 



(A.4) 



An integration by parts gives 
dATiu) 



(3 + aN{u) + {D-l)a I {u - u' + lf~^ Af{u')du'. (A.5) 



du 

where we assumed A^(0) = as there is no sohd present at t = 0. Continuing to differentiate 
gives 



^ ^ - a — 7^ + a{D -1) M{u) 



du^ 



du 
{D-l)\ 



(A.6) 



p-3)!yo 

d^Uiu) d'^Uiu) 



du^ 



a- 
+ a 

+ a 



{u-u' + lf~^J\f{u')du' 
dAf{u) 



dv? 



a{D-l] 



du 





-1)! 


(D 


-3)! 


(D 


-1)! 


(D 


-4)! 



M{u) 



where we used that 



{u-u' + lf-^Af{u')du' 



d r r df(x y) 



(A.7) 



for any function / and constant a. The original integral equation is therefore equivalent to 
a simple D-order differential, 



du^ 



(A.8) 



equation with boundary conditions 

Sf{o) = 0, 

dW{u 



dU{u) 



du 



(A.9) 







du^ 



44 (D- 1)! d^-'Uiu) 



-'^{D dui' 



j = 2...D-l 



The solution is a sum of exponentials, 



D 



ttj exp (A 



.jU), 



(A.IO) 
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where the constants Aj are the roots of 

D 



Using Eg. ( 1A.9p . the coefficients satisfy 



D 



= (A. 12) 

D 

13 = ^AjOj, if D > 1 

D 

a/3 = ^A^^a^-, if D > 2 

i=i 

and so on. In general, solution for the coefficients requires solving this system. For example, 
for the most interesting cases of D = 2 and D = 3 , the solution is 

"-^ ^A.-A,-A. '^^^ ^^-''^ 



Finally, the mass fraction of solid is 

md t'^- 
dt 



s{t) = / m{t,t')^^^dt' (A.14) 



or s{t) = 's{u{t)) with 



siu)= I {u - u' + if ^^^du' (A.15) 

du' 



as GypR^/rriQ = 1. Differentiating gives 

d^, . 



du 



— -^ + D (u-u' + l) — -^du' (A.16) 
du Jq du' 

Substituting from eqn flA.Sp . 

iLj(.) = f5^ + £f^^-^| (A,17) 

du du a \ du 



givmg 



J(u)= (l + -]Af(u)-D^u. (A.18) 
\ a J a 

Insertion of eqn (lA.lOp into eqn (lA.lSp results in eqn (I4.12p . 
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2. Derivation of eqn (I4.17p 



In order to make contact with the phenomenological model, it is useful to develop the 
solution for the mass fraction as a power series. From eqn flA.Gp and the boundary conditions, 
it is easily shown that 



^f{u) = ^ c. 



u 



n=0 



with 



Co = 0, Ci = P, C2 = -a/3, and for > 3 



(A.19) 



(A.20) 



2a 



a 2a 
n n[n — l) n[n — l)\n — 2) 



Cn-3- 



Thus, 



V a I ^-^ a 

^ ' n=l 



(A.21) 



/3m + -/3 (a + 3) + -/3 (a + 3) [a + 2) + ... 
2 6 



We can invert this power series using Lagrange inversion theorem which states that 



u{s) = Y 



n=l 



U 



s{u) 



M=0 



(A.22) 



and results in 



1 (a + 3)^ (2a + 7) (a + 3) ^ 



2/32 6/33 
Finally, to fully eliminate u in favor of s we need 

= /3 + /3 (a + 3) M + J/3 (a + 3) (a + 2) + . . . 
2 

/I / ON ~ 1 (« + 3)^ , 
= /3 + (a + 3)s--^^-^?2 + ... 



(A.23) 



Substitution of this relation into the model, eqn (14. lip , gives the rate-like form, eqn (14.170 . 



3. Size distribution 



Let M[R^ t) be the number of clusters of radius R or smaller. Since the clusters grow 
monotonically, a cluster of radius R at time t was nucleated at some definite time t'{R) < t. 
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Hence, the total number of clusters with radius less than R is the total number of clusters 
minus the number with radius greater than R which is to say the total number minus the 
number already present at time t'{R), 

M {R, t) = [Af (t) - Af {t' (R)) e {R^ (t) - R)] e{R- Ro) (A.24) 

where the step function Q {R — Rq) enforces the condition that there are no clusters smaller 
than Rq and the step function 6 {Rm (t) — R) is required since there is a maximal size 
corresponding to a cluster nucleated at time t = (assuming there are no clusters present 
at t < 0). Hence, when R > Rm, A/{R,t) is simple equal to the total number of clusters 
A/'(t). Now, the problem is to find t'{R). It is sufficient to note that R = R{t,t') = 
(uit)-uit') + l)Ro. Since Afit'iR)) = ^Afiuit'iR))) = ^Af(u + l-^y It then 
follows that 

M 

M{R,t) = Q{R- Rq) — X (A.25) 

U {u{t)) -M(^{t) + l-^^Q {Rm (t) - R) 

Finally, note that Rm is found by taking t' = giving Rm = R(t, 0) = {u{t) + 1) Rq so that 

M 

J\f{R,t) = Q{R- Ro) — X (A.26) 

mo 

J\f {u{t)) - TV" (^M (t) + 1 - -^^ e {(uit) + l)Ro- R) 

This expression simply means that the number of clusters of size R or smaller is the total 
number of clusters created since the time t'{R) at which clusters of size R were created. 



4. No cooperativity 

The limit of no cooperativity, ka = a = 0, is not easy to extract from the general solution. 
A simpler approach is to return to eqn ( lA.Sp which, in this limit, becomes 

^ = Mk^x = Mk^r'^ = ^£P^du 
dt dt kg dt 

or _ 

d^Uit) _ dN{t) _ rriQpR^kndu _ ^du 

dt dt kg dt dt 
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so that Af{u) = [5u. Substituting this into the result for the size distribution gives 

U {R,t) = (3 ( —] Q {R - Rq) (A.28) 
u{t) - U (t) + 1 - ^ ) e {{u{t) + 1) /2o - R) 



Ro, 

The total mass fraction can now be calculated since the mass of a cluster of radius R is 
simply GyR^p and the number of clusters with radius between R and R + dR is ^^j^^^^ dR 
and '^"Ip'*'' = 13—4-0 {R — Rq) 6 {{u{t) + 1) Rq — R). Summing over all radii gives 



M \moJ Jji^ Ro 



Y {{u{t) + 1)^^^ - i) 



D 

In this limit, our result thus reduces to a form that occurs in the theory of crystallization^. 
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